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Abstract 

We investigate both analytically and numerically the evolution of scalar perturbations generated 
in models which exhibit a smooth transition from a contracting to an expanding Friedmann universe. 
If the perturbation equations are formulated as second order equations for either the Bardeen 
potential \& or the curvature perturbation on uniform comoving hypersurfaces £, at best one of 
CN ■ them can stay regular during the transition. We find that the resulting spectral index in the late 

£>v radiation dominated universe depends on which of these two variables passes regularly through the 

transition. The results can be parametrized by the exponent q defining the rate of contraction 
of the universe, or equivalently through the equation of state w = (2 — q)/3q of the background 
fluid. For q ^ — ^ we find that there are no stable cases where both VP and £ are regular during 
the transition. In particular, for < q <C 1, we find that the resulting spectral index is close 
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■ to scale invariant if is regular, whereas it has a steep blue behavior if £ is regular. We also 
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show that as long as q ^ 1, perturbations remain small during contraction in the sense that there 
exists a gauge in which all the metric and matter perturbation variables are small. This work 
has important implications for the current debate concerning the nature of perturbations evolving 
through a collapsing regime into an expanding one: it shows that if in the ekpyrotic model, where 
< q <C 1, the Bardeen potential passes regularly through the transition, this leads to a nearly 
scale invariant spectrum with n = 1 — 2q, whereas in the case of dilaton-driven string cosmology 
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we have the opposite situation. There it is assumed that £ passes regularly through the transition, 
O j' leading to a very blue spectrum of highly suppressed perturbations. 
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I. INTRODUCTION 



Recently, it has been argued that a slowly contracting universe which transits smoothly 
over to an expanding radiation dominated era may lead to a scale invariant spectrum of 
adiabatic density fluctuations [3, Q, 0, 0| . There has been quite an intense debate on the 
question of whether the resulting spectrum after the transition to the radiation era indeed 
is scale invariant or whether it has a steep blue spectrum with n ~ 3. The latter result has 
been advocated mainly by 0, S Q • Gratton et al. have now gone further arguing that the 
ekpyrotic or cyclic scenario is the only robust case where a scale-invariant spectra can be 
found (along with inflation for the case of an expanding universe) 0. 

If the equations of motion governing the transition for the background and the perturba- 
tions were known, this problem could be solved by integrating them numerically. However 
it is likely that the evolution in this high curvature regime will contain full string theory 
and we do not even know whether the variables of the low energy theory are appropriate 
for the description of this regime. One possibility which has been studied in the past is 
the inclusion of first order corrections in the string scale a' and/or the coupling constant g s 
(see, e.g. H 0, 11, El). In this context it has been shown that, within a certain range 



of coefficients for the terms added to the tree level Lagrangian, one can exit from the high 
curvature regime and enter into a radiation dominated phase 0]. In ^3[ the corrections of 
the perturbation equations have been derived and have been solved for dilaton-driven string 
cosmology. Lately, Tsujikawa et al. 0] have used these equations to study the ekpyrotic 
model |4| . Even though they have followed the perturbations through a regularized tran- 
sition, we will argue that their method invariably leads to n ~ 3 and does not allow for 
a decision whether a blue, n ~ 3 or a scale-invariant, n ~ 1 spectrum of perturbations is 
obtained. In this paper we first study a general transition which satisfies relatively mild cri- 
teria and we formulate conditions for the transition which lead to either of the two spectral 
indices. We can decide under which conditions either the spectral index n = 1 as advocated 
by Durrer and Vernizzi 0] or n = 3 put forward in ^| and others is obtained in the ekpy- 
rotic model. This model is a special case of the more general class of contracting universes 
discussed here. 

The rest of this paper is organized as follows: In Sec. |H] we present the aspects of cos- 
mological perturbation theory needed in this paper. We then discuss the transition from 
a contracting to an expanding phase emphasizing the differences between such a transition 
and that associated with the transition between a conventional inflationary phase to radia- 
tion. We also formulate the problem encountered when inferring the perturbation spectrum 
after the transition from the one before the transition. In Sec. IIHI we study the behavior 
of perturbations during a transition and find that the resulting spectral index depends on 
mutually excluding, simple regularity conditions for the transition, which we formulate in 
detail. 

In Sec. IIVI we study numerical toy models for the transition where we exemplify the 
general results obtained in the previous section and analyze the stability of the numerically 
obtained spectral indices. Then we formulate our regularity condition as a theorem and we 
study the amplitude of perturbations. In Sec. IS we comment on the findings in previous 
work and we summarize our results. 
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II. BACKGROUND AND PERTURBATIONS IN SIMPLE CONTRACTING AND 
EXPANDING UNIVERSES 



In this section we repeat the basic equations for the Friedmann background and adiabatic 
first order perturbations. We then discuss initial conditions for perturbations in a contracting 
universe and explain how a different spectral index is obtained depending on the perturbation 
variable used for the transition. 



A. The background 

We consider a Friedmann universe with negligible spatial curvature, which is first con- 
tracting for — oo < 7] < —r] s and hence its (spacetime) curvature is growing. The variable r\ 
denotes conformal time and rj s > is the moment when we enter a high-curvature regime. 
We call this period the "pre-big-bang phase" . At r\ ~ — t] s corrections coming from some un- 
derlying theory, which reduces to general relativity when the curvature is sufficiently small, 
become important and we assume that they regularize the geometry and lead to an expand- 
ing universe at rj > rj s . We also assume that radiation is produced during the transition so 
that for t] > t] s (the "post-big-bang" phase) the universe can be described as a radiation 
dominated Friedmann universe. In the pre- and post-big-bang phases, |?7| ^> t] s , we have 

n 2 = ypa 2 , (1) 

W = -y(p + 3P)a 2 = -^ 2 i±^, (2) 

where the equation of state of the background fluid is w = P/p and n 2 = 8ttG = 2/Mp. 
Mp = l/v47rG = 3.4 x 10 18 GeV is the reduced Planck mass. 7i = a' /a = a = Ha is the 
comoving Hubble parameter. A prime denotes a derivative with respect to conformal time 
rj, whereas a dot denotes a derivative with respect to physical time t, defined by dt = adt]. 
We shall also use 

w ' = 3(w - c 2 )(l + w)H where c\ = ( ] (3) 



V dp 



ad 



q = . (4) 

H l + 3w K J 



is the adiabatic sound speed (if p ^ const, c 2 = P'/p'). We will be especially interested in 
phases during which Pj p = w = const. Then c 2 = const = w and the scale factor evolves 
like a power law, 

1 rj_ 

Vs 

An inflationary phase, defined by 1 + 3w < 0, is thus realized when q < 0. During the 
radiation dominated post-big-bang phase w = 1/3, hence q — 1. 

We will consider a scalar field dominated pre-big-bang phase. During this phase we have 

P=£ + V M , P=£-l%), (5) 



so that 
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For a scalar field, w is constant only for the following possibilities: 

1 if V = , 

w = { -1 if <// = , (7) 

A 2 /3-l if V = V exp(-\ip) , V ^ , 



and, correspondingly, 



1/2 if V = , 

-1 if (f/ = , (8) 
if \/ = K)exp(-A^) , H^O. 



B. Perturbations 



We now discuss perturbation theory in a Friedmann universe (neglecting spatial curva- 
ture, K — 0) with a scalar field or a perfect fluid, like radiation. 

We first consider the linear perturbation equation for the Bardeen potential (see, e.g. [3, 

+ 3W(1 + c^' + [2W + (1 + 3c 2 s )H 2 - TA] $ = . (9) 

This equation is valid for adiabatic perturbations of a fluid, with T = c 2 , or for a simple 
scalar field, with T = 1 (see, e.g. |l6j|). 
If we define the canonical variable 

M P a , , 

u = ; * , 10 



u satisfies the equation [1 
for 



u" + (TAr — 9" /6)u — , (11) 



= a 



or (12) 

o = e 2 = e 1 f%. 

If we restrict ourselves to the case w = c 2 = const, the mass term in Eq. ©, namely, 
2H' + (1 + 3c 2 )H 2 , vanishes by the use of the background Einstein equations (JTJ) and (j2j). 
For these backgrounds which have a oc \r]\ q , where q is given in Eq. (0J), we find 

+ 2H^-^-y' + Tk 2 ^ = . (13) 
Q 

The u equation then simply becomes a Bessel differential equation, 

M » + (T^-^) M = 0. (14) 

The correct normalization to the incoming vacuum at 7] — > — oo determines the initial 
conditions such that (see, e.g. 0|) 



u 
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u = 


rj_ 


-g 


V_ 




Vs 




Vs 



with \i — q + 1/2. Here if/?'* denotes the Hankel function of the second kind and of order 
/i. At "late times" when | A7T7 j <C 1 but still rj <C — rj s we may neglect the k 2 term in Eq. (JUjl 
and find the super-Hubble scale solution 

1+9 

= u A + u B . (16) 

The coefficients A_ and i?_ are determined by the initial solution (|15j) . 

A_ ~ k~ 3/2 (kr] s )- q , B_ ~ k- 3/2 (kr) s ) 1+q , (17) 

hence they have the spectra 

= |A_| 2 fc 3 ~ (/b/ s )" 29 oc k nA -~ x , n A _ = 1 - 2g , (18) 
p B = \B^\ 2 k 3 ~ (£^) 2+2g oc k nB -~ x , n B _ = 3 + 2g . (19) 

Furthermore, ~ |/c?7| 1+2 ' 3 . Hence for g > —1/2, the A mode, u A , dominates at late 

time over the B mode, ub- 

Another perturbation variable often used is the curvature perturbation on uniform co- 
moving hypersurfaces j3| 

H^' + H 2 ^ 

A simple substitution using Eq. © and the background equations yields 

hence on super-Hubble scales, \k/7i\ <C 1, this variable is conserved. For ordinary inflation- 
ary models, it is therefore usually sufficient to compute ( at the time of the Hubble radius 
crossing during inflation to obtain its value in the radiation dominated era. Furthermore, 
since during radiation ( = (3/2)$, this simply gives the Bardeen potential. 
The evolution of ( is closely related to the canonical variable v defined by 



MpaVH 2 - W , 

< ■ (22) 



,r TH 

This variable satisfies the equation 

•" 1 (Tk 2 - z"/z)v = , (23) 



v 



where 



~ _ y — aVH 2 -H' 
Z — Z\ - 



•JfH ' 

or (24) 



z 2 = z 1 J 



dr) 
75 ■ 



Note that the relation between v and ( is v = —MpZ\C,. Equation ()23|) is invariant under 
the "duality" Z\ — > z 2 , in the same way as Eq. (|TT|) is invariant under 9 1 —> 9 2 . 

In [HI it is shown that v appears in the perturbed action as a canonical scalar variable. 
Hence on sub-Hubble scales, \k/H\ ^> 1, it satisfies the initial condition v in = exp(-ikrj) / \fk. 



5 



As before, we now concentrate on the case w = c 2 = const, with the scale factor given in 
Eq. (JU). Then 

Z\ oc a and v = const x M P a( , (25) 

and the constant of order unity depends on q and T. 

During the pre-big-bang phase Eq. (J23)) then also becomes a Bessel differential equation, 

„» + ^_M£_li))„ = o. (26) 

We have set T = 1 throughout the evolution, which should be fine as we are considering 
mainly super-Hubble terms, which satisfy k 2 T <C z"/z during the transition. We have 
confirmed in numerical simulations that allowing T to vary during the evolution does not 
affect the key results we present which concern the spectral index. 
The solution with the correct initial conditions is 



v = V*M H?Hhlh (27) 

with v — 1/2 — q. At "late times" when \krj\ -C 1 but still 77 <C — rj s we may neglect the k 2 
term in Eq. (J26|) and find the super-Hubble scale solution 



v = C_ 


a 


1-? 




a 




Vs 






Vs 



v c + v D . (28) 

The coefficients C_ and D_ are determined by the initial solution (|2~7I) . 

C_ ~ {k^Y^J 2 , D. ~ (krj^t) 1 / 2 . (29) 

The spectra obtained depend on the value of q. One finds 

P c = \C_\ 2 k 3 ~ (kri s ) 4 - 2q ri; 2 oc k nc -~ x , n c _ = 5 - 2q , (30) 

= | 2 fc 3 ~ {kr) s ) 2+2q r)J 2 oc F -- 1 , n D _ = 3 + 2q . (31) 

Here \vc/vd\ — \kr]\ l ~ 2q , hence the D mode dominates for q < 1/2, while the C mode 
dominates for q > 1/2. Finally we want to note that u and v are related via 

Uu/e,)' , (32) 



T 

u = J^ivM • (33) 

It is interesting to note that from the lowest order approximations for u and v given in 
Eqs. (JIEJ), ([2%]) the equivalences (|32~j) and (jHSJ) of w and u cannot be recovered. Only when 
we go to the next order in the term proportional to 9\ or Zi, respectively (or when using the 
full Bessel function solution) does the above equivalence give ua < — ► vc and ub < — > vd 
along with 

n A + 4 = n c , n B = n D . (34) 



6 



C. The problem of a transition from contraction to expansion 



Let us first consider the case where q is in the interval — 1/2 ^ q ^ 1/2. Even though 
q < does not represent a contracting phase, no difference of the following arguments arises 
from letting q decrease until —1/2 (usual inflation has q > — 1). 

Comparing the amplitudes of the modes of u and v, we see that at the transition to the 
expansion phase, \r]\ ~ r] s , we have ua 3> Ub and vp ^> vc for cosmologically interesting 
scales with k ~ l/f]o ~ lCr 30 /?7p. Here t)q denotes the value of conformal time today and 
r/p = Mp 1 . Naively, we therefore expect that immediately after the transition 

u~A^ and v ~ £>_ for - 1/2 < g < 1/2 . (35) 

Since in the radiation dominated era u (x rf and t> oc 77 on super-Hubble scales, we expect 
during the radiation phase 

u ~ A^(n/r] s ) 2 and v ~ D_(r]/r] s ) for kr] <^ I . (36) 

From the relations of m and \1/ as well as w and £ during the radiation dominated phase, this 
gives 

* ~ \[2 A./(r] s M P ) and C ^ £>-/M P • (37) 

For q 7^ —1/2, this naive result is clearly in contradiction with the fact that during the 
radiation dominated era \l/ and ( differ only by a constant since, according to Eqs. (|18J) and 
f!31|) . \1/ would have the spectral index nq, = = 1 — 2q while ( would have the spectral 
index = Ud_ = 3 + 2q. 

For q ^ — 1/2 the B mode of u, ub dominates (for q = —1/2, ua and ub are of the same 
order) and we expect \l/ to have the spectrum = ub_ = 3 + 2q = n£>_ — n^, hence we 
obtain the same spectrum as ( in the radiation era, so that there is no contradiction. 

For q > 1/2 the C mode of v, vq dominates and hence ( actually has the spectrum 
n c, = n C- — 5 — 2g, which is in even worse disagreement with the naively expected spectrum 
for ^ . This contradictory situation is shown on the right hand panel of Fig. ^ Since 
for ordinary inflation q ~ — 1, this problem has never been realized when studying usual 
inflation. 

The simplest possibility which could resolve the issue is to note that the decaying mode 
of u (q > —1/2) during the pre-big-bang phase, ub, has the same spectrum as vd- Hence if 
the w-growing mode during the pre-big-bang phase is entirely converted into the decaying 
mode after the transition and therefore cannot be seen late in the post-big-bang era, we 
expect the spectrum n = 3 + 2q in the radiation era. This argument has been put forward 
in (HI], where the authors have shown that this is exactly what happens if the transition is 
defined by a vanishing jump in the metric and the second fundamental form on the constant 
energy hypersurface. Similar arguments have also been presented in 0,0|. They led these 
authors to the conclusion that the correct spectrum, evaluated sufficiently long after the 
pre — > post transition so that the decaying mode has died away, is n = nq, = = 3 + 2q. If 
this is correct the spectrum of the ekpyrotic model is very blue and in contradiction to the 
observed close-to-scale invariant spectrum. The same argument for the original pre-big-bang 
model of Veneziano , where the scalar field potential vanishes and hence q = 1/2, led 

to the conclusion that the dilaton perturbation spectrum is very blue with n = 4 j2lf . 22 . 



In this argument has been criticized for two main reasons. First, the background 
second fundamental form given byTi/a has to jump, even to change sign, in a transition from 
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contraction to expansion. It then seems quite unnatural to require its perturbation to vanish. 
Second, if the matching conditions are posed on an only slightly different hypersurface, the 
naively expected spectral index, = 1 — 2q is obtained. This "instability" of the index 

= 3 + 2q will also be illustrated in Sec. HVI with numerical studies of a simple toy model. 

The above argument cannot be used if q > 1/2 because u has no mode with spectral 
index 5 — 2q. In this case, agreement can only be achieved if the dominant contribution to 
v, Vr>, is also transferred entirely into the decaying mode so that late after the transition v 
and hence ( still have the spectral index tld_ = 3 + 2q. However, this is not possible: As 
one easily concludes, e.g. from @ or @, a transition on the constant energy hypersurface, 
where the growing mode of u is transferred completely into the decaying mode, preserves (, 
hence ( has the same spectrum after the transition as before, = 5 — 2q which does not 
agree with the spectrum of ^> which in this case is = 3 + 2q. 

Hence if the transition is such that both \1/ and ( correspond after the transition to one 
of their modes before the transition, the obtained spectral index must be n = 3 + 2q. As 
we have shown, this cannot happen for q > 1/2 if the transition is "simple", i.e. does not 
modify the spectrum of either or £. 

If the spectral index after the transition is n = 1 — 2q as promoted in for q > —1/2, 
the variable ( makes a /c-dependent jump at the transition. If n = 5 — 2q is obtained as 
in jl3, ^ makes a fc-dependent jump. Furthermore, if n = 3 + 2q is obtained for q > —1/2, 
the growing mode of \l/ before the transition has to be converted entirely into the decaying 
mode. For q > 1/2 also this no longer helps resolve the problem, and one of the two variables 
C or \I/ must be modified in a fc-dependent way during the transition. 



III. GENERAL SOLUTIONS OF THE PERTURBATION EQUATIONS 
THROUGH THE TRANSITION 

Having explained the problem, but before discussing possible resolutions, let us collect 
some generic facts about a transition from contraction to expansion. Clearly, to have such 
a transition 7i, W and H = (7i/a)' /a = {Ti 1 — 7i 2 )/a 2 have to change sign. Within the 
framework of general relativity (neglecting spatial curvature) this requires p + P < and 
therefore cannot be achieved with a scalar field (with standard kinetic term). If a positive 
spatial curvature is added, the scalar field initial condition can be fine tuned such that 
close to the collapse the curvature term dominates over the scalar field contributions, and a 
transition from contraction to expansion can be achieved with a standard scalar field j23[. 

In this section we want to discuss the problem outlined above without specifying any 
details of the transition. For this, we first discuss the linear second order differential equation 

x" + (k 2 - V x ) x = , V x = — , (38) 
v ' s 

which we have encountered in the previous section. Here the variables (x, s) stand for either 
(u, 6) or (v, z). The factor T in front of the k 2 term is disregarded since it is irrelevant for 
our considerations which mainly concern super-Hubble scales. We notice that Eq. ()38|) is 
invariant under the duality transformation s = Si — » Si j drj/s 2 = s 2 - If s is a power law, 
S\ = I77/77.J 7 , we can set s 2 = Iv/Vsl 1 "" 1 ■ The duality property of Eq. (|3*5|l has been discussed 
in [2J|. If s and 1/s are bounded in the interval [q in ,7]], so that 

00 > C= max fs(r/!) 2 , l/s^x) 2 ) , (39) 
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this equation has the general solution 2(1 Hf 

ctTcos(s, k) + /3Tsin(s, k) 



x 



(40) 



where Tcos and Tsin are defined by 



Tcos(s, k) 



+k 4 



dr] 2 s 2 (r] 2 ) 



Tsin(s,£;) = k I 

J n 



n s 2 (vi) J mn 
s 2 (?7i) 



P 2 dm P 3 

d V 2S 2 (V2) -i^T d VA S 2 ( m )-k 6 --- , (41) 



Vin s2 (Vl) Jr) in 



111 



drj 2 s 2 {r]2) 



'/-> 



dm 



+ k 5 



(42) 



When expressing a given solution in terms of Tcos and Tsin the coefficients a and f3 will 
depend on the initial value t]i n chosen. But as long as s and 1/s are bounded, the sums (JiTj) 
and ()42j) always converge since the terms in this sum are bounded, e.g. by the terms in 
the series expansion for cos [Ck(r] — r] in )] and sin [Ck(r] — r] in )], respectively. Here C is the 
bound from Eq. above. 

To relate this solution with the results of Sec. Ill CI we choose rji n such that k\r] in \ 1, 
but r)i n — r) s . If s obeys a simple power law, s = |?7/?7s| 7 , hence s" / s = 7(7 — l)/r] 2 , we 
obtain to lowest order 

Tcos(s, k) — 1 , 



and 



so that 



x 



Tsin(s, k) = 

(3krj s 



kr]s 



2 7 



1-27 



1-27 



a + 



2 7 



1-27 



1 - 2 7 



1-7 



(43) 
(44) 

(45) 



(For 7 = 1/2 the powers in the Tsin integral becomes a logarithm, but we shall neglect this 
logarithmic correction here.) There are several facts to note at this point: 

• Only two of the three parameters rji n ,a,f] which determine the initial conditions are 
independent. 



As in Eqs. flT2"j) . (J2U), t wo pump fields Si and s 2 yield the same potential V x in Eq. 
If t] in is chosen such that the contributions to the integrals from the lower boundary 
can be neglected, changing s from s\ to s 2 transforms Tcos into fc^Tsin and Tsin into 
/cTcos. 

If s oc rp is a pure power law and again the contributions from the lower boundary can 
be neglected, Tcos oc \fk\t]\ J y {\kr]\) and Tsin oc ^k\r]\Y_„(\kr]\), where v = 7 — 1/2 
and J, Y are Bessel functions. 
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We now define 



Using Eq. (|38j) . we find 



If 7 



(46) 



(47) 



which we can invert to obtain 



-k 2 + (iY s iJ\y' 

Using the latter and Eq. (|38|). the evolution equation for the variable y can be derived 

\ 2 



(48) 



y + 



f--2(i 
s \ s 



y + k 2 y = . 



Let us now introduce 



so that 



r = s 



-r + c 2 s 



ci + c 2 / s 2 d?7 



s 
s 



(49) 



(50) 



(51) 



Equation 



then takes the simple form 

y" + (k 2 - V y ) y = 



(52) 



Note that the y equation obtained from a given x equation depends on our choice of 
s. Since s'[/s\ = s^/s^ but s' x f s\ ^ s' 2 /s2, V y = —s"/s + 2(s'/s) 2 depends on this choice. 
Such a "dual variable" y can also be found if k 2 is modified into T(r])k 2 , the expressions 
just become somewhat more complicated. Choosing x = v and s = Z\, Eqs. (J47|) . ()48|) just 
reproduce the relations (|32|) . ()33|) where y = u and r = 9. During a power law evolution of 
the scale factor, we have Z\ = \i]/i] s \ 9 , z 2 = Iv/Vsl 1 " 9 , @i = \v/Vs\~ q an d ^2 = \v/Vs\ +q ■ 

As we have seen in the previous section, on large scales, \krj\ <C 1, the general solution of 
Eq. ()38)1 is to lowest order of the form 



x = A 


H 


7 


VB 


H 




Vs 






Vs 



1-7 



o{\H 2 ) , 



(53) 



where one obtains from Eq. (J43j) 



1-27 



a = A + B 







B(l - 2 7 ) 
kris 



(54) 
(55) 
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To discuss what might happen during a transition we now assume that for a given s = Si 
or s 2 the solution (|^U|) can be continued through the transition to the radiation era, the only 
effect of the transition being a modification of s which interpolates from 



\v/Vs 



to 



7- 



for 1] <C —Tj s 



(56) 



s = \v/Vs\ 1+ for 7] > rj s , (57) 

without passing through zero. At some time r] s <C rj <C 1/k in the radiation era, we still can 
approximate the Tcos and Tsin integrals by the first term in their series expansion ()41J42j) . 
Integrating the first term in Eq. ()42|) . we obtain 



x 




27- ~ 2 7 h 
l-2 7 + 



Vin 
Vs 



l-27_ 



1-7+ 



Here 



1 - 2 7 _ 

Vs 



-2j~ 



dfj 



(58) 



(59) 



comes from the contribution to the integral during the transition. We always assume that 
the free normalization of s is chosen such that s(—r) s ) = 1. The dimensionless, fc-independent 
constant Ay is then the only quantity that incorporates our ignorance of the true form of 
the transition. Although its typical order of magnitude is O(At) = 0(1), we will mention 
explicit limits on Ay as we go along. With the above expressions for a and (3 we have 



x 



A — B [ A T 



2-2t_ - 274 
1-27+ 



7+ 



B 



1 -2 7 _ 
l-2 7 4 



1-7+ 



(60) 



which is independent of r\ in as it should be. 

In what follows we will study eight different cases and compute the resulting spectra. 
For the first four cases we shall assume that u remains regular throughout. Recalling the 
notation that 9- refers to the collapsing phase and 9 + to the expanding phase in Eqs. fTTf . 
()12)) . we shall consider the following possibilities: = l^/Vsl^ 9 goes over smoothly into 
9\ + = {j]/ri s )~ l (case 1); #i_ goes over smoothly into 9 2 + = (i]/Vs) 2 (case 2); 9 2 - = \r]/f] s \ 1+q 
goes over smoothly into 9\ + (case 3) and 9 2 - goes over smoothly into 9 2 + (case 4). We shall 
then study the equivalent cases for v with 9 replaced by z in Eqs. (|2~3|) . (f2"4"j) . We are mainly 
interested in a contracting pre-big-bang phase, q > 0, but the results derived here are valid 
also for —1/2 < q. 



Case 1: 7 = - 

Here we have A 



-q, 7+ = -i- 

= A. and B = B_ 



hence 



2(2 + q) 



(61) 
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If At is of order unity, or more precisely if 



A T < (k Vs )~^ , (62) 
the resulting spectrum as well as the amplitude does not depend on Ay and we have 

/ . (63 ) 
(if2) 2 (%) 2+2 « ( £J for A; > ^(77) n = 3 + 2g, 



(2fl-2)/(2fl+l) 

fc u (77) ~ r/" 1 ( -1 (64) 



where 

fcu(»?) - ty -1 ( 

is the wave number where we see a kink in the spectrum. For a value of q in the regime 
of our primary interest, — 1/2 < q < 1, the exponent (2g — 2)/(2q + 1) is negative and 
&«(??) <C ?7 _1 especially at late times, r] ^> r) s . Hence, in this case the spectral index relevant 
for the observed anisotropics in the cosmic microwave background (CMB) is n = 3 + 2q, a 
steep blue spectrum. In reaching this conclusion, we have used the fact that in the radiation 
era, 

* - i&" • < 65 > 

/ M \ 2 

|^| 2 A; 3 ~ (^J (£^) 2+2 * iork>k u ( V ), n = 3 + 2q, (66) 

where we have introduced the transition mass scale, M s = 777 1 . For transitions from contrac- 
tion to expansion, q > 0, the amplitude of these fluctuations is therefore far too low to be 
of any relevance for cosmologically interesting scales, k ~ rj^ 1 . (Furthermore, the spectral 
index n = 3 + 2q is not consistent with observations.) 

Case 2: 7_ = —q, 7+ = 2. 

Since 7_ is the same as in case 1, A and B remain unchanged. From Eq. (|6()j). we find 



11 



k~ 3/2 



2 

V 



>, (67) 



so that 

\u\ 2 k 3 ~ (k Vs )- 2 « f^j , (68) 
|*| 2 A; 3 ~ f^j\k Vs )- 2 \ n = \-2q. (69) 

For this result to apply, the condition on A^, Eq. (j62j) must be satisfied. If this case is 
realized and if < q -C 1, a scale invariant spectrum will be obtained. Its amplitude is 
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determined by the transition scale which should be about 5 orders of magnitude below the 
Planck scale. 

Case 3: 7_ = 1 + q, 7+ = — 1. 

According to Eqs. (J53|) and (|TB|) . we now have A = B- and B = A-. This leads to 



u 



k~ 3/2 



(^ s ) 1+ * - (kr) s )- q A T 



2(1 -q) 



-1 



1 + 2q 



(kris)-* 



(70) 



We thus obtain 



\u\ 2 k 



2 1,3 



(kris 



\Wk 



2;„3 



Ms 



-2 Q f V_ 

{kr] s 



-2g 



n = 1 — 2q 



(71) 
(72) 



as in case 2. Here this spectrum is obtained without any condition on A^ having to be 
satisfied, although for the correct amplitude to be obtained, we need A T < (r)/r) s ) 3 /3. 

Case 4: 7_ = 1 + q, 7+ = 2. 

Again, we have we have A = B- and B = A_ and so we obtain 



u 



k~V 2 



(kVsY +g - (kv s r q A 



2(2 + q) 



l + 2q 



(kVs)" 1 ± 



(73) 



with 



\u 



m 2 e 



^ ( ^( At -^)) (a 

V,„,-(a t -^ 



n = 1 — 2q 



(74) 
(75) 



Again, we obtain a scale invariant spectrum if q -C 1, but in this case the amplitude depends 
on the details of the transition given by Ay. 

We now repeat this analysis considering the alternative variable v with 'pump field' s = z\ 
or z 2 . 

Case 1: 7_ — q, 7+ = 1. 

According to Eqs. (|53*j) to (J55J) we have A = D_ and B = C_ , leading to 



^ <^ [(^) g " 1/2 - (^ s ) 1/2 -" (At - 2g)] ( ^ ) + (1 - 2g)(£^) 



(76) 
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MV ~ i " 7 2(% ' )2+2 ' (*) , fM 9 * 1/2 ' (77) 

' iTHkV,)*-* (At - 2<j) 2 (i) for q > 1/2 , 



Id"*- (^) 2 « 2+2 -3 + 2 ?i f^V2, (?8) 

(ft) (^) 4 " 29 (At - 2g) 2 , n = 5 - 2g if g > 1/2 . 

If g < 1/2 we must require A^ < (kr] s ) 2q ~ 1 for our result to apply. Note also that the 
amplitude of the spectrum depends on the details of the transition given by Ay when 
q > 1/2. 

Case 2: 7_ = q, 7+ = 0. 

Again, we have A = D_ and B = C_ , which yields 



fWs { [(^r 1/2 - (kVs) 1/2 - q (A T - 2 + 2g)] - (1 - 2g)(^ s ) 1 / 2 ^ ^ 



(79) 



\v\ 2 k 3 ~ <( ^" 2 (^) 4_2<? (J)' if ? > 1/2 or k > KM , (g0) 
' r/7 2 (A;r/ s ) 2+29 if g < 1/2 and jfe < ^(77) , 

|C| 2 A ,3 (ft)' (^) 4 " 2<? . n = 5 - 2g if g > 1/2 or k > k vi , 

(ft) (^) 2+2? (J)"' , n = 3 + 2g if g ^ 1/2 and fc < ^ . 

For the result to apply when g > 1/2 it requires Ax < T\jr\ s . As in Case 1 of the u field, 
there is a kink in the spectrum with the wave number of the kink for the case g ^ 1/2 being 

which is always smaller than the Hubble radius, r/ _1 , for the relevant values of q, < q < 1/2. 
Only for very small values of q, this kink in the spectrum lies very close to the Hubble radius 
and is not visible. 

Case 3: 7_ = 1 — g, 7+ = 1. 

Here we have A = C_ and B = D_ , hence 

v = v^{((^ s ) 1/2 " 9 -(^ s ) 9 " 1/2 [AT-2 + 2g]) (l\-(kr )s y- 1 / 2 (l-2q)} ,(83) 

H 2 k ^C M ' 2q % ^ >1/2 ' (84) 

V7 2 (kVs) 2+2q (i) [A T -2 + 2g] 2 if g ^ 1/2 , 

| C , V a <(*)>>"■ ' n = 5-2gifg>l/2, 

(ft) (^ S ) 2+2 "[A T - 2 + 2g] 2 , n = 3 + 2g if g ^ 1/2 . 

For g > 1/2 the amplitude of the resulting fluctuations does not depend on the details of 
the transition while it does depend on it for q =Sj 1/2. 
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Case 4: 7_ = 1 — q, 7+ = 0. 

Here again we have A = C_ and B = D_ , so that 



\v\ z k 



2 ;„3 



^|W) 1/2 - 9 - (^ s )^ 1/2 [A r - 2g] + (1 - 2q)(kr )s y- 1 / 2 (jL 
r)- 2 {kr]. 



2 7,3 



ICf* 



where here 



ml 

( Ms.) 

\M P J 



\4-2g 




if 


)2+2g ^ 




2 

if 




-2q 


(*)" 


(kvs) 2 - 






h 

rv V2 


(v) 


~ rf 



n = 5 — 2q if q > 1/2 and fc < fc U2 
n = 3 + 2q if g ^ 1/2 or < k V2 , 



(2-2 9 )/(l-2 9 ) 



For 1/2 < g < 1 and k V2 < r] 1 , a kink from n = 3 + 2g to the final spectrum n 
present in the spectral distribution. 



. (86) 
(87) 

(88) 

(89) 
5 — 2q is 



In Tables HV lIIII we summarize the results of our analysis: 

The spectral index for a transition with regular u and q > —1/2 



case 


7- 


7+ 


kink? 


stable? 


ampl. depends 
on transition? 


n 


1 


-9 


-1 


yes 


no 


no 


3 + 2g 


2 


-9 


2 


no 


yes 


no 


1 - 2q 


3 


1 + 9 


-1 


no 


yes 


no 


1 - 2q 


4 


1 + 9 


2 


no 


yes 


yes 


1 - 2q 



TABLE I: Here we summarize the post-big-bang spectral indices as a function of the pre- and post- 
big-bang exponent of the pump field, if u is regular through the transition. The mild requirement 
on the transition (At) for this to hold is mentioned in the text. 

From Eqs. ()l()j) and ()13j1 it is clear that during the radiation dominated era, inside the 
Hubble radius, kt] 3> 1, the Bardeen potential oscillates and its amplitude decays like 1/rj 2 , 
whereas during the matter dominated era, the Bardeen potential remains constant also inside 
the Hubble radius. Therefore, a change in the spectral index close to the Hubble radius 
crossing is not visible for scales which cross the Hubble radius in the radiation dominated 
era. This remark concerns mainly the kink in case 2 of Table |Hj 

A kink in the spectral distribution arises only in the unlikely situation where the growing 
mode of the pre-big-bang phase is fully converted into the decaying mode after the transition, 
and one has to wait a sufficiently long time [~ i]{k)} for the decaying mode to decay and the 
final growing mode to dominate. It is only in such a situation that the final spectral index 
does not correspond to the naive expectation from the pre-big-bang phase. Finally, we note 
that a kink is always associated with an instability of the spectrum. The issue of stability 
will be discussed in Sec. II VI where we model the regular behavior of u and v through simple 
toy models. There we shall see that a slight modification in the transition can change the 
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The spectral index for a transition with regular v and q ^1/2 



case 


7- 


7+ 


kink? 


stable? 


ampl. depends 
on transition? 


n 


1 




1 


no 


ves 


no 


3 + 2q 


2 


a 





ves 


no 


no 


5-2q 


3 


l-q 


1 


no 


yes 


yes 


3 + 2q 


4 


1-q 





no 


yes 


no 


3 + 2q 



TABLE II: Here we summarizes the post-big-bang spectral indices as a function of the pre- and 
post-big-bang exponent of the pump field, if v is regular through the big bang and q ^ 1/2. The 
logarithmic corrections at q = 1/2 are neglected. The mild requirement on the transition (Ay) for 
this to hold is mentioned in the text. 

The spectral index for a transition with regular v and q > 1/2 



case 


7- 


7+ 


kink? 


stable? 


ampl. depends 
on transition? 


n 


1 


Q 


1 


no 


yes 


yes 


5-2q 


2 


Q 





no 


yes 


no 


5-2q 


3 


l-q 


1 


no 


yes 


no 


5-2q 


4 


l-q 





yes 


no 


no 


3 + 2q 



TABLE III: Here we summarizes the post-big-bang spectral indices as a function of the pre- and 
post-big-bang exponent of the pump field, if v is regular through the big bang and q > 1/2. The 
mild requirement on the transition (At) for this to hold is mentioned in the text. 



spectral index n = 3 + 2q into n = 1 — 2q if u passes through the transition regularly and 
5 — 2q into 3 + 2q if v is regular and q ^ 1/2, correspondingly 3 + 2q into 5 — 2q if v is regular 
and q > 1/2. 

This brings us to one of the key results of this paper, a prediction of the spectral index 
arising from different conditions on q and the regularity of the u and v fields: 



n 



1 — 2q if q > — 1/2 and u is regular and stable, 

3 + 2q if q ^ 1/2 and v is regular and stable, 

5 — 2q if q > 1/2 and v is regular and stable, 

3 + 2q if q < -1/2. 



(90) 



We have not treated explicitly the simple case q < —1/2 above, but this can be done exactly 
along the same lines as the other cases. From Eq. we see that a scale-invariant spectrum 
is obtained for q ~ 1 (standard inflation), or if u is regular and < q <C 1, or if v is regular 
and q = 2. For this latter case, however, we shall see in Sec. IV Bl that perturbations grow 
large during the contracting phase and therefore linear perturbation theory breaks down. 
Furthermore, such a collapsing universe with q = 2 has been shown to be unstable |26j |. 

Clearly, in a transition from contraction to expansion, it cannot be that both u and v 
are regular and stable if q > —1/2. Only in an inflationary transition with q ^ —1/2 do 



1(3 



FIG. 1: Here we illustrate the dominant spectral indices n u (solid) and n v (dashed) on super- 
Hubble scales which are obtained after a "stable" , regular transition of the corresponding variable 
to the radiation era at sufficiently late time. This is also the resulting spectral index for the Bardeen 
potential on cosmologically relevant scales. The left panel shows the indices as a function of 7_ , 
the exponent of the pump field during the pre-big-bang phase. The right panel shows them as a 
function of q, the exponent of the contraction (expansion) law before the big bang. We have used 
Ju- = -Q and 7t,_ = q. 



we find n = 3 + 2q for both u and v. In this case it is expected that both variables transit 
in a regular stable fashion from inflation to the radiation era. The resulting spectral index 
does not depend on the variable with which the calculation is performed. The situation for 
arbitrary values of q is shown in Fig. ^ 



IV. FAST TOY MODEL TRANSITIONS 



In earlier work 0, 0, Q] a transition from contraction to expansion was achieved via 
a combination of first order corrections in the string scale a' and/or the string coupling 
9s- In 03 a modified perturbation equation for v was derived using this framework, and 
a spectral index n = 3 + 2q was obtained. This yields n = 4 for the case dilaton-driven 
string cosmology 0] where q = 1/2, and n = 3 for the ekpyrotic model 3| where q ~ 0. 



However, although calculations have been performed with v, it remains to be shown that u 
cannot pass through the transition regularly (to first order). Unfortunately, even though the 
perturbation equations of 0] are very complicated, they are probably not realistic. It is clear 
that at a time where first order corrections become important, higher order corrections are 
likely to be relevant and the real behavior of the perturbations might differ significantly from 
the results obtained in the work cited above. In this sense, pre-big-bang models including 
first order corrections are only toy models. 

In this section we confirm our generic findings concerning the spectral index associated 
with the relevant u and v fields by numerically solving a simple toy model. We do not 
insist on a physically well motivated transition. Rather we artificially define a regular scale 
factor so that it agrees with a contracting Friedmann universe with contraction exponent q 
at r) <C — r] s and with a radiation dominated universe at rj ^> i] s . In the region in between, 
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-20 20 -2 2 -2 2 




FIG. 2: This figure illustrates the evolution of various background quantities scaled in such a way 
that they are dimensionless (i.e., invariant under a change in rj s ). As functions of ij/r] s , they are 
(1) the scale factor a(rj), (2) the Hubble rate TCr] s and (3) its first time derivative H'%. On the 
bottom line, they are (4) the pre- factor of the comoving wave number T, (5) the potential of the 
perturbation variable u, i.e., V u ?] 2 = (TL 2 — H 1 ) rj 2 and (6) its rescaled square root y|T^J/T rj s . 
The parameters used for these figures are e = 10 -2 and q = 5 • 10~ 2 . 



the scale factor smoothly evolves between contraction and expansion. 
A. The background 

For the exact form of the regularized background scale factor we choose 

«V) = [(v/Vsf + ef 2 , (91) 
q( V ) = qi + (l-t) , (92) 

l( V ) = 1/2 - -arctgfo/^) , (93) 

7T 

where < e < 1. The function l(t)) could be replaced by any C°° function which quickly 
interpolates between 1 for t] —r] s to for rj ^> rj s . Clearly, this universe contracts like \r)\ q 
for r\ — 7] s and expands like a radiation dominated universe a oc rj for t] ^> t] s . Furthermore, 
a, 7i, H! and 7i 2 — H! are all regular, even analytic in the vicinity of the transition, j] = 0. 
The behavior of the relevant background quantities for our model are shown in Fig. El 



B. A regular transition in the perturbation variable u 

We first consider a regular transition for u. In the regime where the scale factor is a 
simple power law, |7/| ^> t] s , the u equation is given by Eq. (fT4"|) . During this regime, the 
u potential is simply V u = q(q + l)/r] 2 . In order to regularize this potential during the 
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transition, we impose V u = 9" /9 throughout, where 



7/2 



0(V) = [(V/Vsf + e] 
l(v) = L l- + (1 - 07+ 



(94) 
(95) 



This ensures us that the pump field 0(7?) remains regular during the whole evolution and 



7+ 



reduces to the power law asymptotic regimes, — ► |7?/?7 S | 7 ~ for 77 <C —r\ s and 9 — > (77/77, 
for 77 ^> 77 s . 

The prefactor of the comoving wave number is T = 1 during the scalar field dominated 



pre-big-bang phase and T 
during the transition via 



1/3 in the radiation dominated era. We regularize T 

(96) 



1 



It is worth mentioning at this point that the results we have obtained appear to be quite 
insensitive to how T(rj) is modelled. To choose the pump field of case 1 of Sec. Ill B[ 9 = 61, 
requires 7_ = — q and 7+ = —1. Similarly, setting 9 = 6*2 corresponds to case 4, with 
7_ = 1 + q and 7+ = 2. To obtain also the cases 2 and 3 we need 9 to interpolate from 9i 
to #2 (case 2) and from # 2 to #1 (case 3), respectively. We can achieve these behaviors using 
our fast interpolating function t(rj) given in Eq. (|9~3j) . The four cases are then obtained by 
the following choices: 



6(1 -q) 


- 1 


(case 1: 


0i- 


-01 , 


-9 - 


> -1) 


t(2 + g) 


+ 2 


(case 2: 


0i- 


-02 , 


-9 - 


> 2) 


t(2 + g) 


- 1 


(case 3: 


02" 


-01 , 


1 +q - 


> -1) 




+ 2 


(case 4: 


02" 


-02 , 


1 + q - 


> 2) 



(97) 



It is easy to verify that the given functional forms have the correct asymptotic behavior. 
Furthermore, they are clearly regular throughout. The numerical results for the w-spectra 
are shown in Figs. El to El The wave number k is given in units of the maximum amplified 
wave number defined by 

k m = max(v/|K|/T) , (98) 
which is of the order of 77" 1 . More precisely we have 



8.9 r/' 1 in case 1 

12.2 in case 2 

3.7 t]' 1 in case 3 

15.1 r]' 1 in case 4 



(99) 
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It is very encouraging that the numerical simulations produce a spectra with precisely 
the shape predicted analytically in the previous section. The spectra are evaluated at 77 = 
Vend — 10 3 ?7 S . As an example, we see that the kink predicted for case 1 is there in the figure 
and is found at the correct position, 

/ n (2fl-2)/(2 9 +l) 

k u = r]- 1 U-J ~ 8 x lCT 9 ??; 1 ~ 10~ 9 £; m . (100) 

One also sees that the spectrum has the correct sub-Hubble radius slope, n su b = n — 4 for 
scales k > 10 — 3 ?7~ 1 which have already entered the Hubble radius at r) en d. 

It is interesting to note that the cases 2 and 3 have roughly the same amplitude while 
the amplitude of case 4, which we expect to depend on the transition, is much higher. 

It is important to investigate the stability of the spectral index n = 3 + 2q of case 1. To 
do this we have slightly modified the potential for this case in the following way: 

1. For r]/r] s G [-0.5,0.5], V U -*V U + KT 1 , [(1) in Fig.Ej. 

2. For r]/r] s G [-0.5,0.5], V U -*V U + 10~ 3 , [(2) in Fig.Ej. 

3. For V / Vs G [-0.1,0.1], V u -> V u + 10~ 3 , [(3) in Fig.Ej. 

The result is shown in Fig El The plain curve represents the original case, V u = Q'[jQ\. 
Adding a tiny constant [~ (9(10 -3 )] to this potential during about a tenth of the duration 
of the transition [curve (3)] already modifies significantly the amplification on super-Hubble 
scales and the final spectral index becomes n = 1 — 2q. Analyzing the growing and decaying 
modes separately, we have seen that, due to the perturbation of the potential, a tiny portion 
of the growing mode during the pre-big-bang phase is converted into the growing mode 
during the radiation era. This is already sufficient for the latter to inherit the naively 
expected spectrum n = 1 — 2q like the other cases. The later we evaluate the spectrum 
the more pronounced becomes the difference from the "pure case 1" spectrum. We expect 
that at very late times, hence very large scales, extremely small differences from the pure 
case 1 potential will have lead to a scale invariant spectrum. We have also tested smooth 
modifications of the potential, like V u — > V u + 10~ 3 exp[— 10 2 (?7/r7 s ) 2 ]. They also lead to the 
same result. The spectra of the cases 2 to 4, however, are stable under small modifications 
of the corresponding potential. 

Finally, in Fig. |S] we show the corresponding spectra for dilaton-driven cosmology. The 
only difference to the previous simulations is that we set q = 1/2 in this model which, in the 
Einstein frame, is a contracting universe with a scalar field with vanishing potential. Again 
we obtain precisely the spectra expected according to the arguments of the previous section. 

C. A regular transition in the perturbation variable v 

In this section we repeat the analysis presented above for the case of a regular equation 
for the variable v. Since the procedure is very close to the one presented above, we can be 
brief here. We again assume that there exists a regular potential V v such that 

v" + {k 2 T - V v )v = . (101) 

In the case of a pure power law scale factor, a = \rj/rj s \ q , we have V v = q(q — l)/r] 2 — z " I z 
where z is either z\ = a or z 2 = a J a~ 2 drj. To regularize the v equation during the transition 
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10 



10' 



10 



10 



10 



P =|u| 2 k 3 



(1) 
(2) 
(3) 




10" 



10"' 



10" 
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10- 2 k/k 10° 10 




FIG. 7: We draw the spectral distribution V u = |u| 2 fc 3 evaluated at rj/r] s = 10 3 (left) and 
r]/r] s = 10 5 (right) for V u of case 1 and the small modifications of the potential detailed in the text. 
Although modifying the potential during the transition may change k m , we kept the same original 
value of k m for all spectra. 



10~ 5 k/ki ° 10" 10 10~ 5 k*10' 
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FIG. 8: We illustrate the spectral distribution Vu 



l k 3 (case i = 1 — 4) as a function of the 



rescaled comoving wave number k/k m The parameters for the simulations are e = 1(T 2 , q = 0.5 
and the spectral distributions are evaluated in the radiation era r)/rj s = 10 3 . 



era we use our interpolation functions t, 7 and T [see Eqs. (jHHJ), (jHSJ) and (jHEJ), respectively] 
and impose 

z(v) = [(v/Vs) 2 + ef 2 . (102) 
To reproduce the four cases for v, we make the following choices: 



7(>7) 



-tq + 1 



1) + 1 (case 1 
iq (case 2 
case 3 
case 4 



z x -»■ zi, 

zi -> z 2 , g -> 0) , 

z 2 -»■ zi, 1 - <? -»■ 1) , 

Z 2 -»• *2, 1 - 9 -»• 0) . 



(103) 
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The resulting spectra are shown in Fig. 




1Cf 10 1CT 5 W° w - w 1CT 5 Wk 10' 



FIG. 9: We illustrate the spectral distribution Vv = |u| 2 fc 3 (case i = 1 — 4) as a function of 
the rescaled comoving wave number k/k m . The parameters for the simulations are e = 10~ 2 , 
q = 5 • 10~ 2 and the spectral distributions are evaluated in the radiation era at r]/r] s = 10 3 . 



Again the spectra are in very good agreement with those obtained with our theoretical 
arguments. This gives us confidence that we really understand what is going on. During the 
radiation era v grows like rj on super-Hubble scales. Inside the Hubble radius the amplitude 
of v remains constant and it begins to oscillate. Neglecting the oscillations, we therefore 
expect for the cases 1, 3 and 4 

P v oc k 2+2q ( Vs / V + k V )- 2 , (104) 

leading to the observed flat spectrum inside the Hubble radius. The spectrum of case 2 
turns from P v oc k^ 2q to P v oc k 2 ~ 2q inside the Hubble radius. For k > 10~ 3 k m the spectrum 
is not very reliable, since it is influenced by the details of the transition. The kink in the 
spectral distribution expected for case 2 is not well visible since q = 5 ■ 1CT 2 is very small 
(see the argument on the previous section). We have repeated this case for a larger value, 
q = 0.3, where the kink can now be seen, as illustrated in Fig. El Finally, note that case 2 
for v is unstable in the same sense as case 1 for u and therefore probably irrelevant. 



V. DISCUSSION 



A. Which variable, u or v? 

We have shown that during a contracting (or inflationary for q < 0) phase where the 
scale factor evolves according to a power law, a oc \r]\ q with q > —1/2, the variables u and 
v acquire a spectral index 

P u oc k~ 2q , n u = 1 - 2q , (105) 

P v oc k 2+2q , n v = 3 + 2q if q < 1/2 , (106) 

P v oc k A ~ 2q , n v = 5 - 2q if q ^ 1/2 . (107) 
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FIG. 10: We illustrate the spectral distribution V v = |v| 2 fe 3 of case 2 as a function of the rescaled 
comoving wave number k/k m . The parameters for the simulations are e = 10~ 2 , q = 0.3 and the 
spectral distribution is evaluated in the radiation era at r)/rj s = 10 3 . We confront our predictions 
for the spectral slopes: above and below the kink at k Vl /k m ~ 10~ 8 — 10~ 7 , the decaying mode is 
still dominant on scales k < k v , . 



We have also shown that, if the corresponding variable transits regularly into the radiation 
era, the spectrum is inherited in this era. But, during the radiation era and on large scales 
(krj «l)n and v are simply related to via [see Eqs. (fTU|) . (|2Uj) and (J22|) ] 

* = (it) = -ir M P ia ~ iv » ( io§ ) 

which therefore has the same spectrum as both u and v. The only possible resolution of 
this contradiction is that either u or v is not regular at the transition and makes a strongly 
fc-dependent jump. But which one? 

We do not know the general answer to this question. It is quite likely to be model 
dependent. Nevertheless, we are able to prove the following statement. 

Theorem. If the perturbed metric remains regular during the transition and if its evolution 
can be described by a second order differential equation for the Bardeen potential then 
a regular ^-variable which satisfies an equation of the form (jllj) can be found. 

Proof. If the metric perturbations are regular, the Bardeen potential which is in general 
given by (see, e.g. 

$ = A + k~ l H(B - k^H'rp) + k~\B - k- l H' T )' , (109) 

is regular too. Here (scalar) metric perturbations are parametrized with the four variables 
A, B, Hl and via 



ds 2 = 

+2H T kjk m \dx > dx r 

and kj = kj/k. 



1 + 2A)dr] 2 - i2Bkjdx j dr) + { (l + 2H L - ^H T ^S 



J"' 



3 

a 2 e~ ik x , (110) 
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During a period governed by general relativity (and in the absence of anisotropic stresses) 
the Bardeen potential satisfies an equation of the form 

vD" + f( v )y> + [T^A; 2 - g(rj)] * = , (111) 

where /, g and T are smooth functions. If the scale factor obeys a simple power law 
evolution, a oc \r]\ q , we have / = 2(1 + q)/r) and g = 0. If / and g remain smooth during the 
transition, we can define 

u = expQ J fdri\v. (112) 
As one easily verifies, this variable satisfies 

u" + (Tk 2 - V u )u = (113) 

with 

V u = g + \f + \f 2 ■ (114) 

V u is well defined, smooth and bounded in any finite interval, so that the differential equation 

9" = VJ (115) 

has two well-defined solutions Q\ and 82 which are the pump fields. Up to an irrelevant 
constant the so-defined variable u coincides with the well known u given in Eq. (JTTj) in the 
asymptotic past r\ <C — r\ s and in the asymptotic future 77 3> r) s . Hence it is our regularized 
variable u. □ 

On the other hand, if ( passes via a regular second order equation through the transition, 
the same theorem leads to a regular f-equation and hence to a spectral index n = 3 + 2q for 
q ^ 1/2. 

This shows again that it is not possible for both \l/ and ( to pass through the transition 
regularly (if — 1/2 < q). This is consistent with the expressions Eq. (120(1 or Eq. (121(1 which 
relate ( and If these equations are also valid during the transition, £ necessarily diverges 
if \1/ is regular and vice versa since 7i and Ti 1 — H? = aH' have to pass through zero in 
a transition from contraction to expansion (see also Fig. |2J). Of course these relations will 
in general be modified during the transition, but according to our results the modifications 
should be such that one of the two variables has to develop a singularity if the other is 
regular. 



B. Amplitude of the perturbations 

During contraction the Bardeen potential grows like ^ oc |^| _ ( 1+2,? ) on super-Hubble 
scales. One actually has 

-(2+4g) 

(kr] s y 2q for q ^ -1/2 . (116) 

Hence ^ may become much larger than 1 for k l/r/ s and \r]\ ~ rj s . Does this imply that 
perturbation theory breaks down during the contraction phase? We show now that this is not 
the case for q ^ 1. First let us note that a quantity relevant to measure the deviation of the 



(M s y 
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\Mp) 


Vs 
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geometry from Friedmann is, for exam ple, the Weyl curvature whose background component 
vanishes. It is well known (see, e.g. [T3|) that the ratio between a typical component of 
the Weyl tensor to a typical component of the background Riemann tensor is given by 
\C/R\ ~ (krj) 2 ty. The geometrical deviation away from Friedmann is thus of the order of 



21,3 



\C/R\ 2 k 



Ms 
M F 



-4? 



;ii7) 



which is always much smaller than 1 on super-Hubble scales for — 1/2 ^ q ^ 1 and |?7| ^r] s . 
Only in a contracting universe with q > 1 do the perturbations become large and hence 
perturbation theory becomes invalid. 

To ensure that perturbations truly remain small, it is necessary to find a gauge in which 
all the metric perturbations are small. We show now that this is so in the off-diagonal gauge 
which also has been used in jilj for dilaton-driven string cosmology. This gauge is defined 
by Ht = Hl = 0. According to Eq. ()109j) . the Bardeen potential is then given by 



q? = A + k'^HB + B') . 
The (ij) Einstein equation implies (see e.g. 

q? = A + k^iHB + B') = —k~ x HB . 
Before the transition, the Bardeen potential is given by 



VH 2 - W 
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(119) 



(120) 



On super-Hubble scales this gives for q > —1/2, using Eq. ()16j) . 



1*1 
with 



M, 



M P 

a = min(2g + 1,2). 



k~ 3/2 (kr) s 



rh 
V 



l+2q 



ii + o(\k V n} 



'121) 



The exponent a > (for q > —1/2) of the first correction to the dominant term can be 
obtained by expanding the Hankel function solution given in Eqs. (llBpiOj) . With Eq. ()119|) 
we then obtain 
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q 



M, 



Mr 



k s \B\ 
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2+2q 
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2q 



{l + 0(\k V \< 



(122) 
(123) 



which remains small on super horizon scales, \kr]\ < 1 during the entire contraction phase, 
rj < -T) s , for q ^ 1 . From A = -k~ l (2HB + B') and Eq. (IT22|l we find that O(A) = 
O^kr/l^B), hence 



k 6 \A\ 



Ms 
Mp 



\krj\ 



2a-2q 



2+2q 



(124) 
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Inserting the value of a given above this becomes 



k?\A\' 



(*); 
(*)' 



\kr)\ 
\krj\ 



2+2q 



4-2q 



2+2q 



2+2q 



for - 1/2 < g < 1/2 
for 1/2 < q . 



(125) 



With B also A remains small on super-Hubble scales as long as q ^ 1 (actually A remains 
small even for 1 < q ^ 2). In this treatment we have neglected the logarithm corrections 
which are present for q = 1/2. Note that it is highly non-trivial that A remains small. This 
is due to the fact that B oc a~ 2 (l + Olkrjl' 7 ) and hence the lowest order contribution to A 
cancels. 

Since the generic form of the perturbed Einstein equations is 0[h+(kr])h+(kr]) 2 h} = 0(A) 
where h and A are typical metric and matter perturbation variables, respectively (see, e.g. 
[TtI| ) . the matter perturbation variables in this gauge are also small on super-Hubble scales. 
More precisely we find from the perturbed Einstein equations in the off-diagonal gauge 



6 



9 -A+*-{k/H)B^O{A) 



2k 



-A + B~B 



m{i + w) 

— In- 1 A' + H~ 2 (2H' + H 2 )A] ~ O(A) . 
3w 



(126) 
(127) 

(128) 



Here 5 and ttl are the density and pressure perturbations, respectively, and i(p + P)k £ is 
the scalar perturbation of the energy flux, T®. To obtain the above results we have used the 
fact that a oc \r]\ q obeys a simple power law and B oc \r]\- 2q as well as A oc |^| max (°' 1 - 2| 3)_ 

Our result that the perturbations remain small hold as long as rj < —r] s , which is the 
epoch when we expect higher order curvature corrections to become important. Typically 
we would expect this to be around the Planck scale. After that, what will happen depends 
on the specific model considered. 



VI. CONCLUSION 

In this work we have analyzed the behavior of scalar perturbations in a transition from a 
contracting to an expanding Friedmann universe. We have shown that, if the perturbation 
equation during the transition can be formulated as second order equations for either \1/ or 
C, regular variables u or v, respectively, can be found. The resulting spectral index in the 
late radiation dominated universe depends on which of these two variables passes regularly, 
and there are no stable cases where both u and v (equivalently \1/ and Q, are regular during 
the transition. 

The resulting spectral index n is given by 

{1 — 2q if \1/ is regular, 
3 + 2q if C is regular and q < 1/2, (129) 
5 — 2q if ( is regular and q > 1/2. 

Our numerical results for the spectral index obtained from a simple toy model are in perfect 
agreement with the more general arguments of Sec. IIII1 
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This result remains valid in an inflationary universe with — 1/2 < q < 0, but has never 
raised any attention since such models cannot produce the observed scale invariant spectrum 
and require w < — 1. For q < —1/2 both variables u and v lead to the same spectral index 
n = 3 + 2q. Therefore, this problem has not been noticed in works on standard inflationary 
models where q ~ — 1. 

We have also shown that, as long as q ^ 1, and we are in a regime where corrections to 
the equations of motion can be ignored, perturbations remain small during contraction in 
the sense that there exists a gauge in which all the metric and matter perturbation variables 
are small. 

We have also argued that the v equation derived from string corrections in ^| has to be 
considered as a toy model, since higher order corrections cannot be neglected in this case. 

Our findings explain that all the literature based on the variable v predicts n = 3 + 2q, see 
especially Refs. 0] and ^H, while when mainly working with u one finds that the spectral 
index n = 3 + 2q is highly unstable and one typically expects n = 1 — 2q. 

This work has the following important implications. 

• // it can be shown that in the ekpyrotic model 0, 0] where < 5 < 1 the Bardeen 
potential passes regularly through the transition, this model leads to a nearly scale 
invariant spectrum with n = 1 — 2q. 

• In dilaton-driven string cosmology we have the opposite situation. There, q — 1/2 
and it has generically been assumed that ( passes regularly through the transition. 
This has been shown to be true to first order in a' in |l3l |. Then the spectral index is 
n = 3 + 2q = 4 leading to a very blue spectrum of highly suppressed perturbation (2l( . 
In this scale invariant spectrum of adiabatic perturbations of the axion field 
can be obtained via the "curvaton mechanism" j^. If, however, \1/ would be regular, 
a red spectrum with n = 1 — 2q is obtained. This would mean a fatal blow for 
dilaton-driven string cosmology, since the perturbations then become very large in 
the radiation dominated era: Since the Bardeen potential is large at the end of the 
pre-big-bang phase and since \C/R\ ~ (kr]) 2 \^\ ~ at Hubble crossing, the Weyl 
tensor becomes larger than the background Riemann tensor at Hubble crossing. Even 
though the Weyl tensor has constant amplitude on super-Hubble scales, the decay of 
the Riemann tensor during expansion leads to a huge increase in the ratio \C/R\. This 
problem only affects red spectra, since for blue spectra \C/R\ at Hubble crossing is 
always smaller than |\I/(A; max )| ~ \C/R\ evaluated at r\ = —r] s for k = /c max - 

Even though we cannot establish from first principles in this work which spectrum dilaton- 
driven string cosmology or the ekpyrotic model have, we nevertheless have formulated suf- 
ficient (but maybe not necessary) conditions on the transition which would allow such a 
decision. 
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